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(54) Measurement of nonlinear properties of an earth formation 



(57) A method of investigating a formation traversed by a borehole includes, measuring at a plurality of 
different static borehole pressures the acoustic Stoneley and/or flexural wave velocities of waves propagating 
through the borehole and formation, and generating an indication of the nonlinearity of the formation by 
processing the velocity measurements. The velocity measurements are processed either by determining a 
fractional change in the measured acoustic velocity and dividing that fractional change by the change in 
borehole pressure to provide frequency dependent acoustoelastic coefficients, or by determining the 
fractional change in the measured acoustic velocity, and subtracting from the fractional change a component 
generated by the borehole fluid and a component due to linear aspects of the formation to provide a nonlinear 
formation component. By processing the velocity measurements at a plurality of frequencies, the nonlinear 
formation components are used to find nonlinear parameters of the formations. The nonlinear parameters are 
then used in conjunction with shear wave velocity information and a database of experimental data, to 
determine the stress in the formation, the strength of the formation, and therefrom, the amount of additional 
stress required to fracture the 
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MEASUREMENT OF NONLINEAR PROPERTIES OF AN EARTO FORMATION 



This invention relates broadly to methods and apparatus for investigating subsurface earth 
formations. More particularly, this invention relates to sonic borehole tools and methods 
for measuring a nonlinear parameter of an earth formation. The invention has particular 
application in using the nonlinear parameter measurements for determining in situ the 
strength of rocks, which information is useful in the production of oil from the formation, 
although the invention is not limited thereto and provides other useful information 
regarding formation parameters. 

The art of sonic well logging for use in determining formation parameters is a well 
established art. Sonic well logs are typically derived from sonic tools suspended in a mud- 
filled borehole by a cable. The tools typically include a sonic source (transmitter) and a 
plurality of receivers which are spaced apart by several inches or feet. Typically, a sonic 
signal is transmitted from the transmitter at one longitudinal end of the tool and received by 
the receivers at the other, and measurements are made every few inches as the tool is drawn 
up the borehole. Depending upon the type of transmitter or source utilized (e.g. dipole, 
monopole), the sonic signal generated by the transmitter travels up the borehole and/or 
enters the formation adjacent the borehole, and the arrival times of one or more of the 
compressional (P-wave), shear (S-wave), Stoneley (tube wave), and flexural wave can be 
detected by the receivers. The receiver responses are typically processed in order to provide 
a time to depth conversion capability for seismic studies as well as for providing the 
determinations of formations parameters such as porosity. 

A method of using a sonic tool in conjunction with the changing of the pressure in the 
borehole is also known. In the defensive publication HI 156 of Siegfried, it is suggested 
that compressional and shear wave speeds, amplitudes, and phase shifts of sonic 
waveforms be measured in the borehole during multiple runs by a sonic logging tool, 
where the borehole pressure is different for each of the runs. According to Siegfried, 

acoustic properties resulting from a change in pressure provides 
an indication of the relative fraciurinc of the formation. 

tar 



While measurements of the compressional and shear waves are useful in quantifying and 
characterizing various parameters of the formation, including fracturing, it will be 
appreciated that to date, there lias been no successful mechanism for making in situ 
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determinations of nonlinear aspects of the formation. For purposes of this invention, it 
should be understood that the term "nonlinear" when used to describe a material relates to 
the fact that a plot of stress versus strain in a material will exhibit some nonlinear behavior. 
Phenomenologically, the strain energy function U(e) of an isotropic elastic solid can be 
written as: 

U(£) = ja,V)e 2 +g(a,p,Y)e 3 0) 
where e is the strain, A and /i are the second order elastic Lame constants, and a, P, and y 
are the third order elastic constants. From equation (1), it will be appreciated that the stress 
(Tis defined by: 

a = dU/de = f{X,n)e + g{a,p,y)e 2 (2) 
where f and g denote general functions of quantities in parentheses. Based on equation (2) 
it is seen that the second order Lame constants are linear constants, while the third order 
constants are nonlinear, and hence measure the nonlinearity of the material. The more 
nonlinear the stress versus strain plot is, the more nonlinear the material is said to be. 
Various manifestations of non-linearity include: the varying of the acoustic velocity in the 
material when the confining pressure changes; the varying of the acoustic velocity in the 
material when the amplitude of the acoustic wave changes; the interaction of two 
monochromatic acoustic beams having different frequencies to create third and fourth 
acoustic beams having the difference frequency and the additive frequency of the two 
incident beams; and evidence of frequencies being generated within the material which were 
not part of any input signal. 

In the oil production industry, rock properties such as sanding, fracturing and borehole 
collapse can be considered to relate to the nonlinear properties of the formation. In each 
case, the strain in the rock catastrophically exceeds that which would be expected from a 
linear stress-strain relationship. As suggested in the parent applications hereto, since the 
less consolidated a formation is. the more nonlinear it is. a measurement of the nonlinearity 
of the formation can provide a measurement of the relative state of the consolidation of the 
formation. As suggested above, whether a layer of a formation is well or poorly 
consolidated, can broadly affect the producibility of the layer and formation, as well as the 
manner in which production is to be carried out. 

It is therefore an object of the invention to provide a measurement of the nonlinearity of a 
formation traversed bv a borehole. 



« 
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According to the present invention, there is provided a method of providing indications of 
formation nonlinearities by using determinations of the acoustic velocity of the formation 
when the borehole traversing the formation is placed under different static pressures. 

It is a further aspect of the invention to provide indications of formation nonlinearities by 
measuring the Stoneley, or flexural wave velocities both before and after a change in the 
static pressure of a borehole, and by subtracting out the nonlinear contribution of the 
borehole fluid from the determinations. 

It is an additional aspect of the invention to determine values of acoustoelastic coefficients 
of a nonlinear formation at a given frequency which are related to the nonlinear constants of 
the formation by measuring Stoneley or flexural wave velocities both before and after a 
change in the static pressure of a borehole. 

Another aspect of the invention is to determine values of nonlinear constants of a formation 
from values of acoustoelastic coefficients determined by measuring Stoneley and/or flexural 
wave velocities both before and after a change in the static pressure of a borehole. 

A further aspect of the invention is to determine in situ a rock strength of a formation by 
utilizing nonlinear constants of the formation which were determined by measuring 
Stoneley and/or flexural wave velocities both before and after a change in the static pressure 
of a borehole. 

An additional aspect of the invention is to utilize determined linear and nonlinear constants 
of the formation to characterize the formation. 

Yet another aspect of the invention is to use determined linear and nonlinear constants 
which characterize the formation to analyze the amount of additional stress which would be 
required to fracture the formation. 

The method of the invention broadly comprises measuring at a plurality of different static 
borehole pressures the acoustic Stoneley and/or flexural wave velocities of waves 
propagating through the borehole and formation, and generating an indication of the 
noniinearity of the formation by processing the velocity measurements. The processing of 
the velocity measurements includes either determining a fractional change in the measured 
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acoustic velocity and dividing the change by the change in the borehole pressure to obtain 
frequency dependent acoustoelastic coefficients which are an indication of nonlinearity, or 
determining the fractional change in the measured acoustic velocity, and subtracting from 
the fractional change a component generated by the borehole fluid and a component due to 
linear aspects of the formation to provide a nonlinear formation component. The nonlinear 
formation components obtained at a plurality of frequencies are used to find nonlinear 
parameters of the formation. In addition, the nonlinear parameters may then be used in 
conjunction with shear wave speed information and a database of experimental data, to 
determine the stress in the formation, the strength of the formation, and therefrom, the 
amount of additional stress required to fracture the formation. 

The system of the invention relates directly to the preferred methods of the invention. The 
system of the invention preferably includes a borehole tool, a borehole pressurizing means, 
and a processing means. The borehole tool includes a monopole and/or dipole source 
which produce dispersive Stoneley and/or flexural waves in the borehole and one or more 
monopole or dipole detectors (e.g. hydrophones) which are spaced from the acoustic 
source and which measure the magnitude of the resulting waves, and means for processing 
and/or sending the measurements uphole, typically via a wireline. 

The borehole pressurizing means preferably comprises a packer type device in conjunction 
with fluid injection means located on the borehole tool for pressurizing a portion of the 
borehole, or a packer type device located on a well head which pressurizes the entire 
borehole. 

The processing means of the invention is typically a VAX computer manufactured by 
Digital Equipment Corporation of Brainard, Massachusetts, or the like coupled to the 
wireline, which receives the measurements provided from the downhole means for 
processing and/or sending and processes the information. In particular, the processing 
means preferably utilizes Prony's method to determine the velocity of the received waves as 
a function of frequency, and then processes that information to determine frequency 
dependent acoustoelast.c coefficients or frequency independent nonlinear parameters of the 
formation. It should be noted that both the acoustoelastic coefficients and nonlinear 
parameters of the formation are functions of in situ stresses in the reference state of the 
formation. The acoustoelastic coefficients are obtained by determining the fractional change 
in the measured acoustic velocities, and dividing the fractional change by the change in 
static borehole pressure. The nonlinear parameters are obtained by determining at a plurality 
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of frequencies the. fractional change in the measured acoustic velocities, and subtracting 
components generated by the borehole fluid and by the linear aspects of the formation from 
the fractional changes in order to provide nonlinear formation components. The nonlinear 
formation components are then processed to find nonlinear parameters (N t and N 2 ) of the 
formation. In turn, the nonlinear parameters may be processed in conjunction with shear 
wave data and a database of experimental data, to determine the stress in the formation, the 
strength of the formation, and therefrom, the amount of additional stress required to 
fracture the formation. 



According to one preferred aspect of the invention, the nonlinear parameters of the 
formation are determined according to an equation: 

( v w„, _ v w,,v) ; = + Q ^ + A v/vl^-H A v/v) /<nMr (3) 

where (v* 0 """ '- v^"*")/ v?™" n AP is the frequency dependent acoustoelastic coefficient, 

C t 

v r ™ v is the velocity of the Stoneley (or flexural) wave in an unpressurized borehole (or in 
a borehole at a given reference pressure), v StoneUy i s the measured velocity of the dispersive 
Stoneley (or flexural) wave at a known pressure, AP is the difference in pressure between 

the reference pressure and the known pressure, Av/v| //nrar is the portion of the fractional 

change in the Stoneley (or flexural) dispersion caused by an increase in the borehole 
pressure that can be calculated from the linear constants of the formation in the ambient 

state, Av/vj^ is the portion of the fractional change in the Stoneley (or flexural) 

dispersion caused by an increase in the borehole pressure that can be calculated from the 
known borehole fluid rionlinearity in the ambient state, C, and C\, are volume integrals 
which are a function of frequency and are calculable in terms of the Stoneley (or flexural) 
wave solution in the ambient state, and N { and N v are nonlinear parameters of the 
formation. 



According to another preferred aspect of the invention, frequencies in the range of 3 kHz to 
6 kHz are utilized in determining the nonlinear parameters. Further, utilizing the nonlinear 
parameters in conjunction with the change in pressure, a multi-frequency inversion of the 
Stoneley or flexural wave velocity dispersions according to AX = B yields a determination 
of two normalized nonlinear constants (c t44 and c, M ) of the formation, where 
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A = \cCcV 



c; ci 



(4) 



X = \N t AP\ 
|/V : AP 



(5) 



B - 



(Av/v| w „ v - Av/v|„ BW - Av/v\ fluid ) fl 



(6) 



and where N { = -c^/c^ and/V, = -c^/c^ , where c M is the formation shear modulus. 

According to yet another preferred aspect of the invention, the stress derivatives of the 
square of the shear velocities in the formation as a function of stress can be determined 
from the determined nonlinear parameters N, and N 3 (or the related nonlinear constants c l44 
and c as ) and from the determinable linear constants of the formation according to: 

P 0 *V,\ __ (2-M)c, 
dS Y ' Y 



66 + i^L + N ^ 



(7) 



Pn dVt _ (vN 2 - N 2 )c„ (N 2 - 2)v c 
dS Y Y 



66 



(8) 



where p 0 is the formation mass density, V n and V u are the two shear wave velocities for 
propagation along direction a first direction ("1") and polarization along perpendicular 
directions thereto ("2" and "3" respectively), S is the uniaxial stress along the one of the 
perpendicular directions, Y is Young's modulus (i.e. a linear parameter of the formation 

given by Y = c 66 (3c i: + c hh )/(c, : + c ft6 ) , and v is Poissons ratio (i.e. another linear 

parameter of the formation) given by v = c I2 /(2c l2 + c 66 ). The linear parameters are 
determined from the compressional and shear wave velocities and from the formation mass 
density. When the slope of the stress derivatives of p 0 V,;, and p^] are positive, the stress 
in the formation is considered to be well below the formation strength, while when the 
slope of the stress derivatives arc negative, the existing stress in the formation is considered 
to be near the formation strength. From experimental data regarding the compressional and 
shear wave velocities of rock formations as a function of uniaxial stress and the shear 
strength, and by utilizing the stress derivatives as calculated, the rock stress is obtained. 



t 
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From the stress in the rock, and from the experimental data regarding the strength of the 
rock, a determination can be made regarding the amount of additional stress which will 
cause the rock to fail (i.e., fracture). 

Additional objects and advantages of the invention will become apparent to those skilled in 
the ait upon reference to the detailed description taken in conjunction with the provided 
figures, in which: 

Figure 1 is a plot of dispersion curve as a function of frequency and velocity for a Stoneley 
wave in a first formation before and after pressurization; 

Figure 2 is a plot of a dispersion curve as a function of frequency and velocity change for a 
Stoneley wave in a first formation before and after pressurization; 

Figure 3 is a plot of dispersion curve as a function of frequency and velocity for a Stoneley 
wave in a second formation before and after pressurization; 

Figure 4 is a plot of a dispersion curve as a function of frequency and velocity change for a 
Stoneley wave in a second formation before and after pressurization; 

Figure 5 is a plot of a dispersion curve as a function of frequency and velocity for a flexural 
wave in the second formation before and after pressurization; 

Figure 6 is a plot of a dispersion curve as a function of frequency and velocity change for a 
flexural wave in the second formation before and after pressurization; 

Figure 7 is a plot showing the frequency sensitivity of coefficients C t and C, to fractional 
changes in the Stoneley wave velocity caused by a borehole pressure increase of unit 
magnitude in the second formation; 

Figure 8 is a plot showing the frequency sensitivity of coefficients D t and D 2 to fractional 
changes in the flexural wave velocity caused by a borehole pressure increase of unit 
magnitude in the second formation; 

Figure 9 is a plot of the plane wave velocities as a function of stress in the second 
formation; 

Figure 10 is a plot of the square of the plane wave velocities as a function of stress in the 
second formation: 

Figure 1 1 is a block diagram of the preferred method of the invention for measuring 
nonlinear properties of an earth formation; and 

Figure 12 is a schematic diagram of the system of the invention which carries out the 
method according to Fig. 1 I. 

Before describing the method and system of the invention in detail, an understanding of the 
theoretical underpinnings of the invention is helpful. In reviewing the theoretical 
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underpinnings, however, various aspects of the invention should be kept in mind. In 
particular, according to the objects of the invention, it is desired to use measurements made 
of the acoustic Stoneley wave or flexural wave velocities in a borehole traversing a 
formation both before and after a change in the static pressure of a borehole in order to 
determine one or more of: acoustoelastic coefficients of the formation; nonlinear parameters 
of the formation; the stress in the formation; and the amount of additional stress which 
would be required to fracture the formation. The propagation of Stoneley and flexural 
waves in a fluid-filled borehole in an ambient condition is adequately described by the 
linear equations of elasticity. However, when the borehole fluid is pressurized, such as 
with the aid of packers at the wellhead, both the fluid and the surrounding formation are 
subjected to biasing stresses. These biasing stresses in the propagating medium result in 
corresponding changes in the effective elastic stiffnesses, mass density and path length 
between two observation (detection) points. Under this situation, the Stoneley and flexural 
wave speeds will change as a function of the increase in the borehole pressure and the 
degree of nonlinearity of the surrounding formation and borehole fluid. The dependence of 
the acoustic wave velocity on biasing stresses in the propagating medium is referred to as 
the acoustoelastic phenomenon, (see Hughes, D.S. and J.L. Kelley, "Second-order Elastic 
Deformation of Solids ', Physics Review 92 pp 1 145-1 149 (1953); Toupin, R.A., and B. 
Bernstein, "Sound Waves in Deformed Perfectly Elastic Materials - Acoustoelastic Effect", 
J. Acoust. Soc. Am. 33, pp. 216-225 (1961); Thurston, R.N. and K. Brugger, "Third- 
Order Elastic Constants and the Velocity of Small Amplitude Elastic Waves in 
Homogeneously Stressed Media", Physics Review 133, pp. A 1604- 1610 (1964)), and 
characteristics of the formation may be described according to acoustoelastic coefficients of 
the formation. 

The acoustoelastic phenomenon may be described by first considering that the propagation 
of a small amplitude wave in an isotropic and homogeneous medium is governed by the 
linear equations of motion. When the propagating medium is prestressed, the propagation 
of such waves are properly described by equations of motion for small dynamic fields 
superimposed on a static bias. A static bias represents any statically deformed state of the 
medium due to an externally applied load or residual stress. The equations of motion for 
small dynamic fields superimposed on a static bias are derived from the rotationally 
invariant equations of nonlinear elasticity by making a Taylor expansion of the quantities 
for the dynamic state about their values in the biasing state. 



♦ 
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When the biasing, state is inhomogeneous, the effective elastic constants are position 
dependent and a direct solution of the boundary value problem is not possible. Under this 
situation, a perturbation procedure can readily treat spatially varying biasing states such as 
those due to a nonuniform radial stress distribution away from the borehole, and the 
corresponding changes in the Stoneley and flexural wave velocities can be calculated as a 

function of frequency. In particular, a modified Piola-Korchhoff stress tensor P in a 

perturbation model yields the following expression for the first-order perturbation in the 
eigen frequency w m for a given wavenumber k,: 

j^v2 a dv - co 2 m j Apvfvyv 

Ao> = * K (9) 

V 

where AP^ = H^vf, (10) 

"«* = + T *P*r + MV< - M») ( 1 1 ) 

^ l ajy0 ~ ~ C qtr0 W 5.6 + C a f yfiAB E AB 

+ W a.L C M + ™ I . M CaMye (12) 
+ W r.P C ajPf} + W p.Q C ajyQ 

T a0 = C ato* W 6. r 03) 
E AB=j{™A.B + »8.a) (14) 

and where the Cartesian tensor notation, the convention that a comma followed by an index 
P denotes differentiation with respect to Xp y and the summation convention for repeated 
tensor indices is used. Although h a/y0 exhibits the usual symmetries of the second-order 

constants of linear elasticity, the effective elastic stiffness tensor H ^ does not have these 
properties as is evident from equation (11). 

Before discussing the various quantities in equations (9)-(14), it should be appreciated that 
the present position of points of material (e.g. formation) may be written as 

X.(&0=£+hL(Xj+iL(X,t) (15) 
where \v denotes the displacement due to the applied static loading of materia] points with 
position vector X in the assumed ambient isotropic state (also called the reference state), 
and u denotes the dynamic displacement vector of material points above and beyond that 

due to the static deformation. The small field Piola-Korchhoff stress P ni in the intermediate 

state can be decomposed into two parts: 
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P = P L + AP ( l6 ) 
where AP is defined by equations (10) - (14), the superscript L denotes the linear 
component of the stress tensor, and the linear component is defined according to: 

P L -c u" 

In equations ( 1 2), (13), and (17), the quantities and are respectively the second 

and third order elastic constants of the material, with Bruggers notation being utilized for 
the nonlinear constants. Brugger, K. "Thermodynamic Definitions of Higher-Order Elastic 
Coefficients", Physical Review 133, No. 6A pp. A161 1-A1612. The various relationships 
among the various notations for the nonlinear constants that appear in the literature are 
known. Kostek, S. et ah, "ThirdOrder Elastic Constants for an Inviscid Fluid", J. Acoust. 
Soc. Am. Vol 94 pp. 3014-3017 (Nov. 1993). In equations (1 1) - (14), T af) , E AB , and w Sy 
denote the biasing stresses, strains, and static displacement gradients respectively. The 
quantity AP denotes the increase in the borehole pressure. 

In equation (9), AP ai are the perturbations in the Piola- Korchhoff stress tensor elements 
from the linear portion P£ for the reference isotropic medium before the application of any 
biasing stresses: p 0 is the mass density of the propagating medium in the reference state; 
and uj represents the eigen-solution for the reference isotropic medium for a selected 
propagating mode. The index "in" refers to a family of normal modes for a borehole in an 
isotropic formation. The frequency perturbations ware added to the elgen-frequency co m for 
various values of the wavenumber along the borehole axis (*,) to obtain the final dispersion 
curves for the biased state. Thus, the modal solution is obtained. Although this method of 
solution is valid for all modes, of particular interest are the Stoneley (m=0), and flexuraJ 
(w= /) modes. 

■ 

An important feature of the perturbation equation (9) is that the volume integral in the 
numerator can be separated into two independent contributions coming from the borehole 
fluid and the formation. The advantage of this feature is that the borehole fluid induced 
velocity change, if known or determined, may be subtracted from the total velocity change 
in order to obtain the acoustoelastic response of the formation. 



* 
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Using the numerical theory set forth above, computational results have been obtained for 
the Stoneley and flexural wave dispersion curves before and after pressuri2ation of a 
fiuidfilled borehole surrounded by an isotropic nonlinear formation. Two distinct 
formations with material properties as listed in Table I below were considered as having a 
borehole of eight inch diameter, with formation I being a "fast" formation, and formation II 
being a moderately "slow" formation: 

TABLE I 







Cm 


P 


v p 






C !I2 




P 


Formation 


GPa 


Gpa 


kg/m 3 


m/s 


m/s 


GPa 


GPa 


GPa 




I 


19.5 


6.50 


2135 


3022 


1745 


-3467 


-1155 


-1541 


-87 


II 


1 1.1 


4.64 


2062 


2320 


1500 


-21217 


-3044 


2361 


-954 



The nonlinearity parameter defined as f3=(3c lJ -hc ut )/2c n is substantially different for the 

two formations. This nonlinearity parameter for rocks varies by several orders of 
magnitude depending on the rock type, porosity, etc., and is substantially larger in granular 
materials as opposed to non- granular materials. 



The borehole fluid is assumed to have a compressional wave speed of 1500 m/s and a mass 
density of 1000 kg/m 3 . The nonlinear behavior of the fluid is usually expressed in terms of 
the parameters A and B appearing in the equation of state, where A is the bulk modulus of 
the fluid. The ratio of parameters B/A for borehole fluid can be taken to be set according to 
B/A = 5, and the corresponding third-order elastic constants are given by c /n = -22.5 GPa, 
c n2 = -13.5 GPa, and c j:j - -9.0 GPa. 



When the borehole pressure is increased by AP, static deformations of the borehole fluid 

and formation are governed by the static equations of equilibrium and continuity of radial 
component of panicle displace formation yields the following biasing displacements, 
stresses, and strains: 

APcr 



.»»>=() (18) 



^ = -^1.^ = ^1.7^=0 (19) 

£rr = _^£1. Ew = ^£1. £z7=0 (20) 
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where a is the borehole radius, R is the radial distance from the borehole axis, and c w is the 
formation shear modulus. Since the stresses and strains in the formation decrease according 
to the square of the radius, formation stresses will have less influence on the Stoneley and 
flexural waves at low frequencies which investigate deeper into the formation; i.e. high 
frequency Stoneley and flexural waves which do not penetrate as far into the formation will 
be more affected by formation stresses which exist due to the borehole. 

It should be appreciated that the volume integrals of equation (9) extend over both the 
borehole fluid and the solid formation. However, the fluid contribution to the first-order 

perturbation in the eigen-frequency co m can be readily computed from the expression: 

J Ac^«;x a dv - col \&p< u ;u;dv 

AQ) ffd =^ t— Vj - (21) 



where 



dA * r i + -Ww. (22) 



Ac f n * = — APS n S<; = 



^ A J + 76 



A p'=^-AP (23) 

and f/ is the mass density of the fluid in the reference state, B is the fluid nonlinear 

parameter appearing in the equation of state (Kostek, S. et al., "Third-order Elastic 
Constants for an Inviscid Fluid", J. Acoust. Soc. Am. Vol 94, pp. 3014-3017 (Nov. 
1993)), and V f denotes the fluid volume. The fractional change in phase velocity resulting 
from changes in the effective stiffnesses and mass density of the material takes the form: 

Av/v = Aa>/G) m (24) 
for a fixed wavenumber. 

Using equation (9) which includes contributions from both the borehole fluid and the solid 
formation, Stoneley wave dispersion curves in formation 1 are shown in Figure 1 both 

before and after borehole pressurization (AP = 13.79 MPa = 2000 psi)- The dashed line in 

Fig. 1 denotes the contribution of the borehole fluid to the total velocity change-as obtained 
from equation (21). It should be appreciated that in the low frequency limit, the change in 
the tube wave speed after pressurization is essentially due to the acoustoelastic effect of the 
borehole fluid. Similarly, in Fig. 2. the fractional change in the Stoneley wave phase 
velocity as a function of frequency is seen for formation 1. As expected. Av/v increases 



- 13 - 



with frequency because the acoustic field becomes more confined to the borehole wall. 
However, the acoustoelastic contribution from the borehole fluid decreases with increasing 
frequency. 

The Stoneley wave dispersion curves and fractional change in the Stoneley wave phase 
velocity for formation II are seen in Figures 3 and 4 respectively. The rock parameters for 
this formation correspond to those of a dry Berea sandstone. Since formation II exhibits a 

higher nonlinearity than formation I, even for a relatively small pressure difference (AP = 

3.4 MPa =500 psi), significant changes are seen in the Stoneley wave speed. The sharp 
increase in the Stoneley wave speed due to a pressure increase at all frequencies for 
formation II results from the high degree of rock nonlinearity, and it will be appreciated that 
the fluid contribution is small compared to the formation contribution. It will also be 
appreciated that the fractional change in the Stoneley wave speed in formation II is 
significantly higher than in formation I, even when the pressure difference is considerably 
smaller. Thus, in formation II, at 4kHz, the Stoneley wave speed of the formation 
increases approximately four percent when the pressure is increased 500 psi . 

The flexural wave dispersion curves and fractional change in the flexural wave phase 
velocity for formation II are seen in Figures 5 and 6 respectively. As is seen in both 
Figures 5 and 6, the fluid contribution is negligibly small due to the relatively low 
nonlinearity of the fluid when compared to the formation. 

The sensitivity of the two formation nonlinear constants N f and N 2 to the Stoneley 
dispersion caused by an increase in the borehole pressure is obtained from equation (9). 
Substituting the linear (second-order) as well as nonlinear (third-order) elastic constants for 
an isotropic medium subject to the biasing stresses caused by an increase in the borehole 
pressure together with the Stoneley wave solution u m (m=0) into equation (9), the change 
of phase velocity at a given frequency can be expressed as set forth above in equation (3): 

( v w,„_ ^.^^ = + CM + Av/v|^ + Av/v|,„ i ,„ (25) 

where C, = \ \J.-—\ (26) 



14 - 



_ c m c \\2 (28) 

~ C 66 

= C " 2 ~ C| " (29) 

4c 66 

and where compressed Voigt's notation has been used for both the linear and nonlinear 
constants. Since there are only three independent nonlinear (third-order) elastic constants 
for an isotropic material, other third-order elastic constants, such as c f44 and c I55 can be 
also expressed in terms of c nn c //2 , and c m such as: 

c, 55 = (l/4)(c nr c ii2 ) (30) 
c i44 = (\/2)(c u3 -c m ) (31) 
It can be seen from equations (28) -(31) that N, = -c^c^ „ and N 2 = -c /55 /c 66 . It can also 
be seen from equations (26) and (27) that C, and C, are frequency dependent, and their 

frequency dependence for a unit increase in borehole pressure (i.e., AP = -IN/m ) is shown 

in graphical form in Fig. 7. 

It is clear from Fig. 7 that the magnitude of C 2 , is approximately zero over a frequency 
band of 0 to 2 kHz. Therefore, inversion for the nonlinear constant A r 2 of the formation 
from the fractional changes in the Stoneley wave velocities in this frequency band would 
result in large errors in its estimated values. Thus, a preferred frequency band for the 
inversion of N t and N : is one where both the values of C, and C, are significantly large. 
However, it is known in the art that a conventional sonic tool with sources and receivers 
placed on the borehole axis is not very efficient in recording Stoneley wave signals at 
frequencies much higher than about 10 kHz (see, Kurkjian, A.L., "Numerical Computation 
of Individual Far-Field Arrivals Excited by an Acoustic Source in a Borehole", Geophysics 

• ■ * 

Vol. 50, No. 5. pp. 852-866 (May 1985)). Using this information as well as the sensitivity 
analysis of Fig. 7, an optimal frequency band for processing the Stoneley wave velocity 
changes caused by a borehole pressurization for the determination of N l and N 2 will 
preferably be between approximately 2.5 kHz and 7 kHz and more preferably in a band 
between 3 kHz and approximately 6 kHz. 

In solving equation (25) for A' and A' : . it will be appreciated that C, , C : , A\7 v ;W and 

AW vl. are all determinable, and that the term (v*"" 1 "" - v% mH ")/v%~ u * which relates to 

the acoustoelastic coefficients of the formation is determinable. In particular, since the 
Stonelev wave solution docs run have any azimuthal dependence. J r J : . J, and J v of 
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equations (26) and (27) can be expressed in terms of line integrals as set forth in Appendix 
A hereto. Using the line integrals which utilize the shear and compressionaJ wave velocities 
of the formation, C, and C 2 are found. Also, the portion of the fractional change in the 
Stoneley dispersion caused by an increase in the borehole pressure can be calculated from 
the linear constants of the formation in the ambient state according to: 

Av/V \unear = J A I 2(0^ J N (32) 

whero the integral J 4 is expressed as a sum of two terms J 4l and J 42 which are likewise 
defined in Appendix A hereto. Further, the portion of the fractional change in the Stoneley 
dispersion caused by an increase in the borehole pressure that can be calculated from the 
known borehole fluid nonlinearity in the ambient state is given by: 

Av/v U" =y s /y * (33) 
where J 5 is expressed as an integral as set forth in Appendix A hereto. Thus, by measuring 
the fractional change in the Stoneley velocity ((v w " v - v%~ Uy )/ v 5 r ^\ (which relates to 
the acoustoelastic coefficients) the only unknowns in equation (25) are N t and N 2 . 
However, a solution for /V, and N 2 is possible by utilizing equation (25) at a plurality of 
discrete frequencies and conducting a multi-frequency inversion as described in more detail 
hereinafter. 

The sensitivity of the two formation nonlinear constants N t and N 2 to the flexural 
dispersion caused by an increase in the borehole pressure is likewise obtained from 
equation (9). Substituting the linear (second-order) as well as nonlinear (third-order) elastic 
constants for an isotropic medium subject to the biasing stresses caused by an increase in 
the borehole pressure together with the flexural wave solution u m (m=I) into equation (9), 
the change of phase velocity at a given frequency can be expressed as: 

( v f«-- _ w *- ) , v , r « = Dt „ i+ + Ay / V |^ +Av/ v|/ _ (34) 



Where D. = \ — f /, - 1 

" 2co-J x { - 2) 



(35) 



(36) 



N * = •'•> (37) 



•V, = 



(38) 
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In solving equation (34) for N. and JV : , it will be appreciated that D } , D 2 , Av/ v\ fluid and 

Av/ v|, lW are all determinable, and that the term (v flmral - v£" rai )/ v% urai which relates to 

the acoustoelastic coefficients of the formation is determinable. In particular, it will be 
appreciated that since the flexural wave solution exhibits azimuthal dependence, /,, l 2 , I 3 
and I N can be expressed in terms of surface integrals as set forth in Appendix A hereto. In 
addition, the portion of the fractional change in the flexural dispersion caused by an 
increase in the borehole pressure that can be calculated from the linear constants of the 
formation in the ambient state is given by Av/ v\ Unear = IJI N where the integral l 4 is 
expressed by: 

/ 4 = fl 2 AP(/ 4l + / 43 ) (39) 

where I 4J and 1 42 are defined in Appendix A hereto. Similarly, the portion of the fractional 
change in the flexural dispersion caused by an increase in the borehole pressure that can be 
calculated from the known borehole fluid nonlinearity in the ambient state is given by 

& vlv \fiuid = hi 1 ** where the integral I s is defined in Appendix A hereto. Thus, by 
measuring the fractional change in the flexural velocity ((v* w - v^)l v£ w ) (which 

relates to the acoustoelastic coefficients) the only unknowns in equation (34) are N t and N 2 . 
However, as was the case with the Stoneley dispersion, the solution for jV, and N 2 is 
possible by utilizing equation (34) at a plurality of discrete frequencies and conducting a 
multi-frequency inversion as described in more detail hereinafter. 

It will be appreciated from equations (35) and (36) that D, and D 2 are frequency dependent, 

and their frequency dependence for a unit increase in borehole pressure (i.e., AP = -lN/m) 

is shown in graphical form in Fig. 8. Here it is clear that the magnitude of D, and D 2 are 
essentially zero from 0 to 3 kHz. This implies that the inversion for the nonlinear constants 
N { and N 2 of the formation from the fractional changes-in the flexural wave velocities in 
this frequency band would result in substantially large errors for their estimated values. A 
preferred frequency band for the inversion of jV, and N : is one where both the values of D, 
and D, are significantly laree. However, it is known in the an that a conventional sonic tool 
with dipole sources and receivers placed on the borehole axis is not very efficient in 
recording flexural wave sicnals at frequencies either below about 1 kHz or much higher 
than about 10 kHz. iSce Sinha. B.K. et ah. "Borehole Flexural Modes in Anisotropic 
Formations-. Geophysics. Vol. 59. No. 7. pp. 1037-1052 (July 1994)). With this 
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information, as well as the sensitivity analysis of Fig. 8, an optimal frequency band for 
processing the flexural wave velocity changes caused by a borehole pressurization 
approximately 6 kHz, and more preferably between 4.0 kHz and 5.75 kHz. 



As suggested above, estimation of the formation nonlinear constants may be carried out 
from a multi-frequency inversion of the Stoneley and/or flexural wave velocity dispersions 
caused by a borehole pressure increase. Assuming that fractional changes in the Stoneley 
wave velocities at two frequencies /, and f 2 are available for a borehole pressure increase 
aP, the inversion process is formulated as AX = B according to equations (4) - 6) set forth 
above: 



A=\c{ l cV\ 



k /2 C: /2 | 



(40) 



N.AP 



(41) 



B = 



< Av/v| w/fv - Av/v\ iinear - &vfv\ fluid ) fx 
(^/v| w/fv -Av/v| Wflr - Av/v fiuid ) n 



(42) 



It should be appreciated that although equations (40) - (42) relate to a two frequency 
inversion, three or more frequencies can be utilized in the inversion. 

Table II below contains the input data of the Stoneley velocity differences (AP = 500 psi) at 

four frequencies which are grouped in pairs, and the nonlinear constant parameters of the 
second formation are calculated according to equations (25) and expanded equations (40) - 
(42). The actual values of the formation nonlinear constant parameters are shown in the 
parentheses. 



TABLE n 

Multi-frequency inversion of Stoneley Dispersion 



t kHz 




B 




N : 


1.6S8 


0.0168 


0.0149 


582 


977.5 



2.180 


0.0219 


0.0201 


(582) 


(979) 


2.537 


0.0255 


0.0237 


582 


979 


3.574 


0.0347 


0.0331 







It is seen that the values of the nonlinear parameters N { and N 2 calculated from the above 
equations based on obtained velocity information for both pairs of frequency data agree 
closely with actual values for the formation. 



As with the Stoneley dispersion, changes in the flexural dispersions caused by borehole 
pressurization can be subjected to a multi-frequency inversion to estimate the nonlinear 
constants of the same formation. However, the flexural dispersion data is preferably taken 
at a moderately high frequency band (approximately 3 to 4 KHz), as the low frequency 
flexural dispersion data exhibit negligibly small acoustoelastic effects. Thus, Table III 
below contains the input data of the flexural velocity differences at four frequencies which 
are grouped in pairs, and the nonlinear constant parameters of the second formation are 
calculated according to equation (34) and expanded equations (40) - (42), with the matrix 
of equation (40) being slightly modified to account for "D" values instead of M C" values. 

TABLE m 



Multi-frequency Inversion of Flexural Dispersion 



f kHz 


AW v 


B 


N, 


N 3 


3.412 


0.0155 


0.0153 


582 


980 


3.995 


0.0257 


0.0253 


(582) 


(979) 


4.054 


0.0266 


0.0263 

* 


586.6 


978.6 


1.404 


0.0322 


0.0318 







Again, it is seen that the calculated values for the nonlinear parameters N t and /V, for both 
pairs of frequency data agree closely with the actual values for the formation. 

While nonlinear parameters A', and N : by themselves provide additional useful information 
regarding the formation, the parameters can also be used to help determine the type and 
strength of the rock formation which is being investigated, as well as the quantitative stress 
in the formation. According to theory, the plane wave velocities for waves propagating 
along the X, direction in an isotropic medium subject to homogeneous biasing normal 
stresses (T) and strains (E) arc obtained from the equations of motion for small dynamic 
fields superimposed on a bias: 



» 
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PoVZ =c„ +T U +(2c„ +c lll )£ 11 + c 1I2 (£^ + E 33 ) (43) 
PoV,i = ^ 6 + T u + c,„£, , + (2c w + ClJJ )E 22 + c 144 £ 33 (4 4) 
A>K> = + 7J, + c !55 £„ + c l44 £ 22 + (2c 66 + c 15s )£ 33 ( 45) 
where p 0 is the mass density of the medium, V„ denotes the plane wave velocity in the 
reference state for propagation along the X, direction and polarization along the X 3 
direciion. If the specimen is in the form of a rod with the uniaxial stress applied along the 
rod axis (i.e., stresses normal to the rod axis are assumed to be zero), stress derivatives of 
shear waves propagating normal to the rod axis and polarized parallel and normal to the 
stress direction can be approximated by: 

as y y (46) 

PoW (vN 2 -N 2 )c m (/V, - 2)^ 

as y + ? (47) 

where N x = -c 144 / C66 andJV, = -c l5J /c M , N, and N 2 are the normalized nonlinear constants 
of the formation, and v and Kare the respectively Poisson's ratio and Young's modulus of 
the formation in the reference ambient state. The Poisson's ratio and Young's modulus are 
a function of the second order constants of the formation and can be expressed by: 



v = 



C I2 



2 ( C !2 +C 66) 



(48) 



v_ c 6a(3c,, +2c 66 ) 



(50) 



On the other hand, if the specimen is long along the propagation direction, a plane strain 
approximation normal to the propagation direction X, is an appropriate assumption. In this 
case, when the specimen is subjected to a uniaxial stress S along the X 2 direction, stress 
derivatives of shear waves polarized parallel and normal to the stress direction are given by: 

as { -' y 1 r 

Therefore, stress derivatives ofp.,1^ and p 0 V x \ can be readily approximated from either 
equations (46) and (47. or ,50) and (51). and typically, the stress derivatives for a rock 
sample will fall between these approximations. The stress derivatives are functions of the 
formation nonlinear constants c !4J and r,„ .via nonlinear parameters N, and N,), and the 
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linear constants v and Y, and provide two additional formation parameters which can be 
used in addition to parameters such as density, and p-wave and shear wave velocities to 
identify the lithology of the formation. As set forth above, the nonlinear constants can be 
estimated from the inversion of the changes in the Stoneley or flexural dispersions over an 
appropriate frequency band caused by a known change in the borehole pressure. 

It should be appreciated that the stress derivatives set forth in equations (46), (47), (50) and 
(51) may be viewed as the slopes of curves which relate the square of the plane wave 
velocities V , 2 and V u to the uniaxial stress S. Thus, if it is possible to obtain different 
formation samples, subject them to different known confining pressures, stress them with 
known stresses, and measure and plot the shear (plane) wave velocities as a function of 
stress, a database of curves of different formations can be obtained where the stress 
derivatives set forth above will relate to the slope of the plotted curves. 

With that in mind, experimental data was obtained for the compressionai and two shear 
wave velocities in uniaxially stressed cylindrical rock samples of the second formation 
(parameters of the second formation being set forth in Table I above). The rock samples 
were subjected to a uniaxial stress along the cylindrical axis (with no confining pressure), 
and the compressionai and shear waves (the shear waves being polarized parallel and 
normal to the stress axis) were excited and detected by transducers mounted normal to the 
stress axis, with measurements being carried out with increasing stress until the material 
fractured. The results of the experiment are seen in Figures 9 and 10, with Figure 9 being a 
plot of the plane wave velocities as a function of stress, and Figure 1 0 being a plot of the 
square of the plane wave velocities as a function of stress. 

Figure 10 reveals that when the stress derivative (i.e. curve slope) is positive, the existing 
stress in the material is significantly less than the yield (i.e. failure) stress; whereas when 
the stress derivative is negative, the material is closer to failure. In fact, as the negative 
stress derivative increases in magnitude, the rock becomes closer to failure. 

Substituting the estimated values of the formation nonlinear constants N, and N : obtained 
from the multi-frequency inversion of the Stoneley and flexural dispersions before and after 
borehole pressurization and as shown in Tables II and III, the stress derivatives of p 0 V; : 
and p,,Vp for the assumed reference state given by 5=0 were calculated from equations 
(46) and (47). When these stress derivatives are plotted in terms of the quantities shown in 
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Fig. 10, the resulting slopes are given by lines labeled A, and B, respectively in Fig. 10. On 
the other hand, when equations (50) and (5 1 ) are employed in the calculations of the stress 
derivatives, the resulting slopes are given by lines marked A, and B 2 . The agreement 
between the calculated slopes is good in view of the different plane stress and plane strain 
assumptions as well as the experimental uncertainties in the data. 

With the theory as set forth above, the preferred method of the invention for measuring 
nonlinear properties of an earth formation is shown in Figure 1 1 in a block diagram format. 
At step 100, with the borehole at a first known pressure (typically ambient), one or both of 
Stoneley and flexural waves, and preferably both shear and compressional waves generated 
by the borehole tool are measured by receivers or detectors of the borehole at a first location 
in the borehole. The various velocities of the different waves are calculated at step 102 
utilizing known methods. For example, the Stoneley wave and flexural wave velocities 
(i.e., the velocity dispersions as a function of frequency) are calculated utilizing Prony's 
method, while the compressional wave velocity may be obtained from the detected 
compressional headwave. Then, at step 104, the borehole is pressurized to a second known 
pressure. At steps 106 and 108, the Stoneley and flexural waves, as well as shear and 
compressional waves are again generated and measured by the borehole tool at the first 
location, and the velocities (and velocity dispersions) of the different waves are calculated 
at the new pressure. Utilizing the Stoneley (or flexural) wave velocity dispersions, at 1 10, 
the acoustoelastic coefficients are calculated as a function of frequency according to 

V V ) lv r.f *P (or the equivalent flexural wave equation). The frequency 

dependent acoustoelastic coefficients so obtained can be used to provide a qualitative 
indication of the nonlinearity of the formation. The acoustoelastic coefficients can also be 
used as described below in determining the nonlinear constants N, and N 2 of the formation. 

With the Stoneley wave and/or flexural wave dispersion curves, an indication of the 
nonlinear constants N, and N, of the formation are obtained at step 1 12 utilizing equations 
(25) (27). (30) - (33), and <40) - (42). In particular, with a knowledge of the borehoie 
radius, the change in pressure, the linear constants of the formation (typically as determined 
from the shear and compressional wave velocities), and the mass density and nonlinear 
parameters of the fluid, and using equations (26), (27), (30) - (33). (40) - (42) and the 
equations in the Appendix hereio. equation (25) is solved for both N, and /V,. Substeps 
involved in solving equation (25) include, for each of at least two frequencies (preferably 
located in the optimal frequency hands discussed above with reference to Figs. 7 and 8). 
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determining a fractional change in the measured acoustic velocity at step 1 12a, and at step 
1 12b subtracting from the fractional change a component generated by the borehole fluid 
and a component due to linear aspects of the formation to provide a frequency dependent 
nonlinear formation component (B). Then, utilizing an inversion process AX = B at step 
1 12c according to equations (40) - (42), values are obtained for the nonlinear parameters N, 
and N 2 . It should be appreciated that in lieu of utilizing the inversion process, equation (25) 
can be solved at two or more different frequencies to yield two or more values for 
frequency dependent nonlinear formation components which arc each equal to 
C,N,+C 2 N @ . By substituting determinable values for C, and C 2 as a function of frequency, 
the equations can be solved simultaneously for the two unknowns N, and N 2 which are 
frequency independent. Regardless, once obtained, parameters N, and N 7 may then be used 
at step 1 14 to find the third-order elastic constants c tS5 and c, 44 according to equations 
N t = ~q 44/^66 and/*/, = -c^/c^ . 

If both sides of equation (25) are divided by the change in pressure (AP), the left side of 
equation (25) evolves into the acoustoelastic coefficients of the formation as discussed 
above with reference to step 1 10. Thus, if desired, the acoustoelastic coefficients of the 
formation can actually be utilized in the inversion process or in a system of simultaneous 
equations to find the indication of the nonlinear constants N, and N 2 and two of the three 
third-order elastic constants therefrom. 

The nonlinear constants of the formation N, and N : provide an excellent indication of the 
relative consolidation of the formation, as well as formation strength. It is also known that 
a crossplot of rock linear elastic constants expressed in terms of the compressional and 
shear velocities together with the mass density helps in identifying the rock lithology or 
mineralogy (see e.g., Pickett. G.R.. "Acoustic Character Logs and Their Application in 
Formation Evaluation". Journal of Petroleum Technology, Vol. 15 p.650 (1963). This 
crossplot. also known as Pickett's crossplot helps in the identification of rock lithology 
(e.c, identification of rock type as limestone, dolomite, or clean sandstone, etc.) in terms 
of a constant ratio of compressional wave velocity to shear wave velocity of different 
masnitudes. However, in more poorly consolidated rocks, the data tend to diverge from a 
constant ratio resulting in a difficult interpretation of rock lithology. The nonlinear 
constants N, and A', add another dimension to such crossplots which serves to differentiate 
rock material properties on a finer scale. In particular, nonlinear elastic parameters are much 
more sensitive 10 the chan-jc.s in the material structure and initial stress than conventionally 
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measured linear elastic moduli. (See, e.g. Yu, lrina et al., "Comparison of Linear and 
Nonlinear Elastic Moduli for Rocks by Use of a Granular Medium Model", Abstract for 
127th Meeting of Acoustical Society of America, Cambridge, Massachusetts June 6-10, 
1994). 

According to the invention, the determination of the nonlinear constants N f and N 2 of the 
formation can be further utilized in conjunction with a database of experimental data in 
order to determine the stress in the formation, the strength of the formation, and therefrom, 
the amount of additional stress required to fracture the formation. In particular, at step 120, 
a preliminary determination of formation lithology is made. The determination can be made 
from previous knowledge regarding the formation obtained from formation cores, or by 
neutron type logging tools, gamma ray type logging tools, etc. A determination of 
formation lithology may also be made utilizing determinations of density, shear wave 
speeds, compressional wave speed, and other information (including the nonlinear constant 
values N t and N 2 ). (See, e.g. Castagna, J.P. et al., "Relationships Between 
Compressional-Wave and Shear-Wave Velocities in Elastic Silicate Rocks", Geophysics 
Vol. 50, p. 571 (1984); Han, D.H., "Effects of Porosity and Clay Content on Wave 
Velocities in Sandstones", Geophysics, Vol. 51, p. 2093 (1986); Kowallis, B. et al., 
"Velocity-Porosity-Clay Content: Systematics of Poorly Consolidated Sandstones", J. 
Geophysical Research Vol,89, p. 10355 (1984). Then, at step 122, the nonlinear constants 
N t and N 2 are used in conjunction with a knowledge of the linear constants of the formation 
in accord with equations (46) - (51) in order to provide indications of the stress derivatives 
of p 0 V^ and p 0 ^ 2 v At step 124. the values of the stress derivatives of p Q V* 2 and p 0 VJ; 
(i.e. the slopes) are used to define tangents to database plots of the square of the plane 
wave velocities as a function of stress (such as in Fig. 10) of known formations of the 
same or similar lithology at the appropriate confining pressure of the formation. Preferably, 
the database plots include plots which relate the square of the shear velocities to stress in 
different rock samples (lithologies) at different confining pressures, and the slope value 
information is applied to an appropriate database plot. Based on the particular plots and the 
values of the indications of the stress derivatives, the actual stress in the formation is 
determined at step 126. as the slope of the curve defines where along the curve the 
formation stress is located. In addition, because each database plot is obtained by taking a 
formation sample and stressing it until it fails, the strength of the formation is known. By 
subtracting at step 12S the actual stress in the formation from the strength of the formation, 
a determination is made as to the amount of additional stress which would be required to 
fracture the formation. 
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Different techniques for matching the value of the stress derivatives to the slope of the 
database curves (plots) at step 124 can be utilized. A preferred technique is to take each 
curve, select, e.g., four local data points of the curve, fit a cubic polynomial to the four 
local data points, and take the derivative of the polynomial to obtain a slope at that location 
of the curve. The slope and an associated formation stress value are then stored. This 
process is repeated for additional points of the same curve (overlapping the points if 
desired), until the curve is represented by a plurality of different slopes with associated 
different formation stress values. For each lithology, different curves are obtained to 
represent the situation at different confining pressures. Thus, in matching the stress 
derivatives to the slope of the database curves at step 124, a determination of the borehole 
depth is first made in order to determine the confining pressure in the formation. Based on 
that confining pressure and a knowledge of the lithology, the appropriate curve represented 
by the plurality of different slopes with associated formation stress values is chosen from 
the database. Then, the determined in situ stress derivative value is matched to one or more 
of the stored slope values for the appropriate curve and the actual stress in the formation is 
taken as the stress associated with that stored slope value(s) or interpolated therefrom at 
step 126. 

It will be appreciated by those skilled in the art, that one or more logs of values of the 
acoustoelastic coefficients, the nonlinear parameters such as JV, and /V 2 , the rock strength, 
the stress, and the additional stress required to fracture the formation can be made by 
conducting the method of the invention at different locations in the borehole. 

The preferred system 200 for carrying out the method of the invention is seen in Fig. 12. 
The system 200 of the invention preferably includes a borehole tool 220, a borehole 
pressurizing means 250, and a processing means 290. The borehole tool 220 is seen in 
Fig. 12 suspended in a borehole 201 by means of a wireline 203 and a winch 205 as is well 
known in the art. The preferred borehole tool 220 is a tool such as the DSI tool of 
Schlumberger (DSI being a trademark of Schlumberger) which includes a plurality of 
acoustic detectors 222, and one or both of a monopole source 224 and a dipole source 226. 
The monopole source 224 provides a Stoneley wave as well as a compression^! headwave 
in all formations, and also provides a shear headwave in fast formations. The dipole source 
226, on the other hand provides a flexural wave, and the shear wave arrival time can be 
found as the low frequency limit of the flexural wave dispersion arrivals. Thus, in accord 
with the invention, it is possible in fast formations to use only a monopole source, as the 
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Stoneley wave can be used to find the acoustoelastic coefficients, while the Stoneiey wave 
information In conjunction with the shear wave and compressionaJ waves information can 
be used to find the nonlinear parameters N t and N 2 , as well as the actual stress in the 
formation and the amount of additional stress which would be required to fracture the 
formation. On the other hand, in slow formations, use of both the monopole source and 
dipole source is preferred as the iow frequency limit of the flexural wave dispersion curve 
can be used to provide the shear wave speed. Where both monopole sources 224 and 
dipole sources 226 are used, they are preferably pulsed at different times. In addition, 
where both monopole 224 and dipole sources 226 are used, the acoustic detectors 222 
should include both monopole and dipole detectors, or detectors such as four hydrophones 
placed ninety degrees apart from each other and capable of working in both dipole and 
monopole modes. In the dipole detector mode, signals from diametrically opposite 
hydrophones of the four hydrophone arrangement are subtracted to yield the pressure 
gradient at the borehole axis associated with borehole flexural waves, while in the 
monopole detector mode, the signals from the diametrically opposite hydrophones are 
added to yield the average pressure at the borehole axis associated with Stoneley wave 
propagation along the borehole axis. 

The borehole tool 220 also preferably includes a downhole processing means 230, packer 
means 240, and fluid injection means 245. The downhole processing means typically takes 
the form of a microprocessor and associdted circuitry which is coupled to the detectors 222 
of the borehole tool 220. The downhole processing means obtains voltage information 
relating to pressures seen by the detector, and processes the information for transmission 
uphole via the wireline 203. If desired, the downhole processing means can determine the 
shear and compressionaJ wave velocities, as well as processing the flexural or Stoneley 
wave information via Prony s method. Typically, however, the heavy mathematical 
processing of Prony s method is accomplished uphole by the uphole processing means 290 
such as a computer which is coupled to the wireline 203. 

The packer means 240 of borehole tool 220 preferably comprises first and second inflatable 
packers located on opposite ends of the borehole tool 220 such that the sources and 
detectors of the borehole 220 are located between the packers. The packers are preferably 
arranged so that they can seal with the borehole wall and thereby seal off a portion of the 
borehole. The fluid injection means 245 includes a fluid storage means 247, and pump 
means 249 for injecting the fluid into the borehole. When the packers are inflated and in 
sealed relationship with the borehole wall, the borehole (and formation) between the 
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packers can be pressurized by injecting fluid into the borehole. With the borehole being 
pressurized, the monopole and/or dipole sources are energized so that data is obtained for 
processing in accordance with the method invention. 

It will be appreciated that the uphole processing means 290 processes the information 
received from the downhole tool 220 via the wireline 203 according to the method of the 
invention (with steps 102, and 108 - 128 of Fig. 11 providing an extremely high level 
flowchart of the programming of the uphole processing means 290). In particular, utilizing 
Prony's method, the velocities of the received Stoneley or flexural waves are determined by 
the processing means as a function of frequency for the borehole in at least two different 
states of pressurization (unpressurized being one of the possible pressurization states). The 
compressional wave velocity and shear wave velocity are also determined according to 
known techniques. The processing means 290 then processes the Stoneley or flexural wave 
information to determine acoustoelastic coefficients according to the method invention set 
forth above. Also, using either the acoustoelastic coefficients or the fractional changes in 
the measured acoustic velocities at a plurality of frequencies, in conjunction with other 
information as described with reference to Fig. 11, the processing means 290 provides 
nonlinear parameters of the formation. Further, using database information and the 
nonlinear parameter information, and as described in Fig. 11, the processing means 290 
can determine the stress in the formation, the strength of the formation, and therefrom, the 
amount of additional stress required to fracture the formation. 

There have been described and illustrated herein systems and methods for measuring 
nonlinear properties of an earth formation utilizing a sonic borehole tool. While particular 
embodiments have been described, it is not intended" that the invention be limited thereto, as 
it is intended that the invention be as broad in scope as the art will allow and that the 
specification be read likewise. Thus, while a particular borehole tool embodiment was 
described with packers and fluid injection means, it will be appreciated that other means for 
pressurizing the borehole could be utilized. For example, rather than providing the borehole 
tool with packers and fluid injection means, a packer type device can be located on a well 
head in order to pressurize the entire borehole. Also, while the borehole tool was described 
as a wireline type tool, ii will be appreciated that the invention can be embodied in 
measurement-while-drilling (MWD) type tools. Further, while specific equations were 
utilized in providing a theoretical basis for the invention, it will be appreciated that the 
invention is not tied exactly thereto, and that other acoustoelastic coefficients and nonlinear 
parameters relating to aspects of the formation can be obtained by the method of the 
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invention. Likewise, while the invention was described primarily in terms of utilizing a 
Stoneley wave with the equations directed primarily to Stoneley waves, it will be 
appreciated that flexural waves can be utilized to obtain the same results. It will therefore be 
appreciated by those skilled in the an that yet other modifications could be made to the 
provided invention without deviating from its spirit and scope as so claimed. 

APPENDIX A 
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where the changes in the fluid bulk modulus AA, and mass density p f caused by an 
increase in the fluid pressure are calculated from the expressions 
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where and p, are the fluid and formation mass densities, respectively; and w/) 

denote the flexural wave solution in the borehole fluid, whereas ( u r , and u. ) are the 
corresponding solution in the formation. 
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where the changes in the fluid bulk modulus AA, and mass density p x caused by an 
increase in the fluid pressure are calculated from the expressions 
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CLAIMS 



A method of measuring a nonlinear property of an earth formation surrounding a 
borehole using a sonic logging tool comprising a monopole or dipole source and a 
sonic detector, the method comprising: 

a) using the sonic logging tool to generate and measure a first frequency dependent 
velocity of a Stoneley wave or flexural wave at a location in the borehole at a first 
pressure; 

b) pressurizing the location in the borehole to a second pressure; 

c) using the sonic logging tool to generate and measure a second frequency 
dependent velocity of a Stoneley wave or flexural wave at the location in the 
borehole at the second pressure; and 

d) determining the nonlinear property of the formation from the first and second 
frequency dependent velocities (v) according to (v - v rtf )lv rtf AP, where AP is the 
difference between the first and second pressures. 

A method as claimed in claim 1 . comprising determining a frequency dependent 
acoustoelastic coefficient in order to determine an indication of the nonlinear 
property of the formation. 

A method as claimed in claim 1 or 2, comprising determining the fractional change 
in the wave velocities and subtracting from the fractional change components related 
to the borehole fluid and to linear aspects of the formation to provide a nonlinear 

i 

formation component. 

A method as claimed in claim 3, wherein the nolinear formation component is 
provided substantially according to the equations: 

L.siomirs _ ^wfr, j / y ^«r„> = c, A/, + C 2 N 2 + Av I + Av I v|, M 



.tear 



or 



\iear 



i . ii.-.u,*i _ ,..'^-~')/i>' ,u "" = D t N t + D Z N 2 + Avlv\ l1uut + Ai7v|,„ 
wherein Avlv\,. is the component related to the linear aspects of the formation, 

is the component related to the borehole fluid. N, and A\ arc nonlinear 
constants of the formation. C. and C, are volume integrals which are a function of 
frequency and are claculable in terms of a Stoneley wave solution in an abicnt state 
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of the formation, and £>, and D 7 are volume integrals which are a function of 
frequency and are claculable in terms of a flexural wave solution in an abient state 
of the formation. 

5 A method as claimed in claim 4, further comprising measuring a first shear wave 
velocity and a first compressional wave velocity at the first pressure and using these 
velocities to determine C n and measuring a second shear wave velocity and a 
second compressional wave velocity at the second pressure and using these 
velocities to determine C? 

6 A method as claimed in claim 4 or 5, comprising solving the equation for N, and N 2 
by obtaining values for v s """*\ v*°""'\ C, and C 2 for a first frequency {/,) and a 
second frequency (f 2 ). 

7 A method as claimed in claim 6, wherein the first or second frequency is chosen 
from a frequency band of 2.5 kHz to 7 kHz. 

8 A method as claimed in claim 5, 6 or 7, wherein solving the equation comprises 
conducting a multifrequency inversion according to AX = B where 

A = \c/ , cr\ 
\cf 2 a n \ 



X = \N l AP\ 
|/V,AP| 



B = 



< Al '/v| W) „-Ai//v| /(n<or - Av/v\ fliiid ) fl 
(Av/v| JnMI ,,, > -Ai-/v|,,, irar - Av/v flud ) f2 



A method as claimed in claim 4. comprising solving the equation for N, and N, by 
obtaining values for ,."<"•"". y*""* 1 . D, and D, for a first frequency (/,) and a 
second frequency </,). 
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10 A method as claimed in claim 9, further comprising measuring a first shear wave 
velocity and a first compressional wave velocity at the first pressure and using these 
velocities to determine £>,, and measuring a second shear wave velocity and a 
second compressional wave velocity at the second pressure and using these 
velocities to determine D v 

1 1 A method as claimed in claim 9, wherein the first or second frequency is chosen 
from a frequency band of 3.5 kHz to 6 kHz. 

12 A method as claimed in claim 9, 10 or 11, wherein solving the equation comprises 
conducting a multifrequency inversion according to AX = B where 

A=\D( X D{ X \ 
D{ 2 Dl 2 \ 

X=\N { AP\ 
\N 2 AP\ 

B= (Av/v ^w-Av/v^-Av/vl^),, 
(Av/v ^w-Av/v^ - Av/v^)^ 

13 A method as claimed in any of claims 4 to 12, further comprising obtaining a 
database of experimental data for a plurality of other formations relating stress to 
squares of shear wave velocities, and utilizing the nonlinear constants and the 
database of experimental data to find the stress in the formation and/or the additional 
stress required to fracture the formation. 

14 A method as claimed in claim 13, further comprising choosing a stress-shear 
velocity squared curve from the database relating to a formation of substantially 
similar lithology to the formation surrounding the borehole, and finding a value 
relating to the slope of the curve in order to find a point along the curve by utilizing 
the nonlinear constants. Young's modulus for the formation in an ambient reference 
state. Poisson\s ration for the formation in an ambient reference state, and values of 
linear parameters of the formation, the point along the curve providing an indication 
of stress in the formation. 
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A method as claimed in claim 14, wherein the value relating to the slope of the 
curve is found according to one of the equations 

Pogjii A 2 ~ g^jfg , {N x + Aggg 

as Y Y 

Po av,] = (vn 2 - n 2 ) C66 (yy^Ks 

as Y Y 

or 

Y Y 

where V l2 and are shear wave velocities in the formation having lithology 
substantially similar to the formation under investigation along a first direction ("I") 
and polaraized along perpendicular directions thereto ("2" and "3" respectively), S 
is the uniaxial stress along one of the perpendicular directions, Y is Young's 
modulus, and v is Poisson's ratio. 

Apparatus for measuring a nonlinear property of an earth formation surrounding a 
borehole, comprising, a sonic logging tool comprising a monopole or dipole source 
and a sonic detector, and borehole pressurizing means for changing the pressure in 
the borehole at the location from a first pressure to a second pressure; wherein the 
tool is used to generate and measure first and second second frequency dependent 
velocities of a Stoneley wave or flexural wave at the location in the borehole at the 
First and second pressures, processing means being provided for determining the 
nonlinear property of the formation from the first and second frequency dependent 

velocities (v) according to (v - v ref )/v ref AP, where AP is the difference between the 
first and second pressures. 

Apparatus as claimed in claim 16, wherein the borehole pressurizing means 
comprises first and second inflatable packers on the borehole tool with the source 
and receiver located therebetween, the packers being inflatable to contact and seal 
with the borehole wall. 
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18 Apparatus as claimed in claim 17, wherein the borehole pressurizing means 
includes fluid storage means for storing fluid in the tool, and fluid injection 
means coupled to the fluid storage means for injecting the fluid into the 
borehole when the first and second packers are inflated . 

Apparatus as claimed in claim 16, 17 or 18, wherein the tool generates and 
measures a First shear wave and a first compressional wave at the first 
pressure, and a second shear wave and a second compressional wave at the 
second pressure, the processing means using the measurements of the first and 
second shear waves and compressional waves to determine the nonlinear 
property of the formation. 

Apparatus as claimed in any of claims 16 to 19, further comprising a database 
of experimental data for a plurality of other formations relating stress to 
squares of shear wave velocities, the processing means utilizing nonlinear 
constants and the database of experimental data to find the stress in the 
formation and/or the additional stress required to fracture the formation. 

21 Apparatus as claimed in claim 20, wherein the processing means chooses a 
stress-shear velocity squared curve from the database relating to a formation 
of substantially similar lithology to the formation surrounding the borehole, 
and finds a value relating to the slope of the curve in order to find a point 
along the curve by utilizing the nonlinear constants, Young 1 s modulus for the 
formation in an ambient reference state, Poisson's ratio for the formation in 
an ambient reference state, and values of linear parameters of the formation, 
the point along the curve providing an indication of stress in the formation. 
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